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We report quantum and semi-classical calculations of spin current and spin-transfer torque in a 
free-electron Stoner model for systems where the magnetization varies continuously in one dimension. 
Analytic results are obtained for an infinite spin spiral and numerical results are obtained for realistic 
domain wall profiles. The adiabatic limit describes conduction electron spins that follow the sum 
of the exchange field and an effective, velocity-dependent field produced by the gradient of the 
magnetization in the wall. Non-adiabatic effects arise for short domain walls but their magnitude 
decreases exponentially as the wall width increases. Our results cast doubt on the existence of a 
recently proposed non-adiabatic contribution to the spin-transfer torque due to spin fiip scattering. 



I. INTRODUCTION 

The use of electric current to control magnetization 
in nanometer-sized structures is a major theme in the 
maturing field of spintronics.^ An outstanding exam- 
ple is the theoretical prediction of current-induced mag- 
netization precession and switching in single domain 
multilayers^ and its subsequent experimental confirma- 
tion in spin-valve nanopillarsn^ The physics of this spin- 
transfer effect is that a single domain ferromagnet feels 
a torque because it absorbs the component of an inci- 
dent spin current that is polarized transverse to its mag- 
netization. The same idea generalizes to systems with 
continuously non-uniform magnetizationi^ This realiza- 
tion has generated a flurry of experimental^ and theo- 
retical work^'SiSiifliii focused on current-driven motion of 
domain walls in magnetic thin films. 

The experiments cited just above employ Neel-type do- 
main walls with widths w « 100 nm. This length is very 
large compared to the characteristic length scales of the 
processes that determine the local torque. ^^'^^ Therefore, 
it is appropriate to adopt an adiabatic approximation 
where the spin current is assumed to lie parallel to the 
local magnetization^^'- Surprisingly, the adiabatic pre- 
diction for the current dependence of the domain wall 
velocitySiiiii^ agrees very poorly with experiment. This 
has led theorists£*2i2ii2iii to consider non-adiabatic ef- 
fects and experimentersiSiii to study systems with do- 
main wall widths that are much shorter (w « 10 nm) 
than those studied previously. 

Two groups^i^ have studied the effect on domain wall 
motion of a distributed spin-transfer torque represented 
by a sum of gradients of the local magnetization with 
constant coefficients. For a one-dimensional magnetiza- 
tion M(a;), the torque function can be written in terms 
of two vectors perpendicular to the magnetization 

Nst(a;) = cid^M -\- caM x d^rM. (1) 

In general, the coefficients ci and C2 are functions of posi- 
tion. The well-established adiabatic piece of the torque is 



the first term in Eq. ^ with a constant coefficient. Con- 
sistent with usage in the literature, we call all deviations 
from the adiabatic torque non-adiabatic. Any contribu- 
tions of the second term are then called non-adiabatic. 
Zhang and LiS derive a contribution along this second 
direction in Eq. JQ) from a consideration of magnetiza- 
tion relaxation due to spin-flip scattering in the context 
of an s-d exchange model of a ferromagnet. ^'^•^^ Their 
arguments lead them to the estimate C2/C1 « 0.01. The 
authors of Ref. 10 report that a similar value of C2 /ci pro- 
duces agreement with experiment when Eq. |^ is used 
in micromagnetic simulations. 

In this paper, we study the applicability of Eq. Q to 
a free-electron Stoner model with one-dimensional mag- 
netization distributions of the form 

M(a;) — Ms (sin 9 cos 0, sin 9 sin (f), cos 6) . (2) 

In Eq. (O, the polar angle 6{x) is measured from the pos- 
itive z axis, 0(2;) is the azimuthal angle in the x-y plane, 
and Ms is the saturation magnetization. We begin with 
the spin current and spin torque for an infinite spin spiral 
with constant pitch. This system turns out to be per- 
fectly adiabatic; the torque is described by Eq. ^ with 
C2 = 0. The same is true for realistic domain walls of the 
sort usually encountered in experiment. Non-adiabatic 
effects appear only for very narrow walls. In that case, 
the torque is non-local and cannot be written in the form 
Eq. 1^. The non-adiabatic torque decreases exponen- 
tially as the wall width increases for all realistic domain 
wall profiles. Finally, our analysis casts doubt on the ex- 
istence of a non-adiabatic contribution to the torque due 
to spin-flip scattering proposed recently by Zhang and 
Li.9. 

The remainder of the paper is organized as follows. 
Section describes our Stoner model and the methods 
we use to calculate the spin current and spin-transfer 
torque. Section fllll reports our results for an infinite spin 
spiral and Section [IVI does the same for one-dimensional 
domain walls that connect two regions of uniform magne- 
tization. Section Ivl relates these calculations to previous 



work by others. Section rVTI discusses the effects of scat- 
tering. We summarize our results and offer some con- 
clusions in Section rVIII Two appendices provides some 
technical details omitted in the main body of the paper. 



II. MODEL & METHODS 

The free electron Stoner model provides a first approx- 
imation to the electronic structure of an itinerant ferro- 
magnet. The Hamiltonian is 



The adiabatic approximation^ to the spin dynamics leads 
to a spin current density that is proportional to the local 
magnetization, Qad(2;) oc M(a;). This means that 



Nad(a;) cx d^M{x) 



(8) 



A main goal of this paper is to study the extent to which 
the spin-transfer torque associated with real domain wall 
configurations satisfies Eq. 0. 



A. Quantum 
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where cr = {ax, ay^cTz) is a vector composed of the three 
Pauli matrices and fi — ghe/2mc. The magnetic field 
Box(a;) is everywhere parallel to M(a;) but has a constant 
mag nitudeiiS, That magnitude is chosen so the Zeeman 
splitting between the majority and minority spins bands 
reproduces the quantum mechanical exchange energy in 
the limit of uniform magnetization: 



E,^ = 2^|B,,| = h^kl/m. 



(4) 



If Ep = h'^kp/2m is the Fermi energy, the constant ks 
in Eq. fixes the Fermi wave vectors for up and down 



spins, kp and kp, from 



kp 



(5) 



Given Bcx(a;), we use both quantum mechanics and a 
semi-classical approximation to calculate the spin accu- 
mulation, spin current density, and spin-transfer torque. 
The building blocks are the single-particle spin density 
s±(x,kx) and the single-particle spin current density^^ 
Cl±{x,kx) for an up/down (±) electron with longitudi- 
nal wave vector kx- 

Summing over all electrons and using the relaxation 
time approximation gives the non-equilibrium majority 
and minority spin density s-t (x) and spin current density 
Q±{x) in the presence of an electric field Ex: 



s±{x) 



/±(k-^*)-/±(k) 



/±(k-^*)-/±(k) 



s±(a;, kx)^^^ 
(6) 

Q±{x,kx)(fk 



Our use of the function /±(k) — Q{kp — |k|) implies 
that the distribution of electrons outside the region of 
inhomogeneous magnetization are characteristic of the 
zcro-tcmperature bulk.^" We shall expand this point and 
comment on the general correctness of © in Sec. I VII 

The sum s{x) = s+{x) + s_(a;) is the total spin ac- 
cumulation (spin density) and Q{x) = Q+(a;) + Q-(a;) 
is the total spin current density. Finally, the distributed 
spin transfer torque is^^ 



In light of Eq. the exchange magnetic field that 
the enters the Hamiltonian in Eq. (PJ is 

Bex (a;) = Bex (sin cos 0, sin 6* sin 0, cos 6') . (9) 

For a given energy, the eigenfunctions for a conduction 
electron with wave- vector k take the form 5'±(r,k) = 
tp±{x,kx)e'^^^y'^^'^'^ where the spinor tlj±{x,kx) satisfies 



dx"^ 



COS u 
e^"^ sin( 



In this expression. 



(10) 
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and ± refers to majority/minority band electrons. The 
single-electron spin density and spin current density arei^ 

s±{x, kx) = V'±,a(a;, kx)cra.,f3i'±,fj{xi kx) (12) 



and 



Cl±{x,kx) 



■— ^ Im['0J^„(a;, kx)(Ta,p—'tp±,i3{x, kx)]. 



N,t{x) = -dxQix) 



(7) 



(13) 

As a check, we used this formalism to calculate the 
equilibrium (zero applied current) spin density Scq(a;) and 
equilibrium spin current density Qcq(a;) for a magnetiza- 
tion distribution chosen arbitrarily except for the con- 
straint that |M(a;)| be uniform. The densities Scq(a;) and 
Qeq(a;) are obtained by retaining only the second term 
in square brackets in Eq. 0.^^° The lines in Fig. QJa) 
are the Cartesian components of the imposed M(a;). 
The solid dots in Fig. ^a) show that the spin density 
Scq(a;) is parallel to M.{x), as expected. Similarly, the 
electron-mediated spin-transfer torque should equal the 
phenomenological exchange torque density discussed by 
Browni2i This is confirmed by Fig. ^b), which shows 
that Neq(x) is proportional to M x M"(a;) for the M{x) 
shown in Fig. ^a). 



B. Semi-classical 

A semi-classical approach to calculating the spin cur- 
rent density is useful for building physical intuition. 
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FIG. 1: Equilibrium (zero-current) results: (a) Cartesian 
components of an arbitrarily chosen magnetization M(a::) 
(lines); Cartesian components of the calculated spin density 
Seq(x) (solid dots); (b) Exchange torque (solid line) and cal- 
culated spin-transfer torque (solid dots). 
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FIG. 2: (a) Cartesian components of an imposed magnetiza- 
tion M(a;) used in the other panels; (b) comparison of quan- 
tum (lines) to semi-classical (solid dots) calculations for the 
Cartesian components of the spin density s(a;); (c) same com- 
parison for the Cartesian components of the spin current den- 
sity Ql{x). 



Accordingly, we write an equation of motion for the 
spin density of every electron that contributes to the 
current. This idea has been used in the past, both 
semi-quantitatively'* and qualitatively^ Our derivation is 
based on the behavior of an electron with energy E that 
moves along the x-axis through a uniform magnetic field 



i?cxZ. The wave function for such an electron is 



b e**-'-^ 



where 
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2mE/h^±k%. 



(14) 



(15) 



We compute the spin density s{x,E) and the spin current 
density Cl{x^E) for this electron using the right sides of 
Eq. H12|) and Eq. H13|l . respectively, with %l]± — > ip. 

It is straightforward to check that the components of 
these densities transverse to the magnetic field satisfy the 
semi-classical relations 

Qx = Sx{v) and Qy = Sy{v), (16) 

where (w) is the velocity 

^^^-^^ = ^^'^- 

Moreover, the transverse components of the spin density 
satisfy 



dx 
da; 



(18) 



These equations are the components of the vector equa- 
tion 



ds ^ 

dx (fc) ° 



(19) 



where Box = z. 

We now make the ansatz that all three Cartesian com- 
ponents of the semi-classical majority and minority spin 
densities s±{x,kx) satisfy Eq. (|19|l when the direction 
of the magnetic field varies in space. Specifically, if 



Box(a;) 



^{x), we suppose that 



dx (k) ' 



k,) X Bex(x), (20) 



where k+ and fc_ for s+{x,kx) are defined by Eq. H15|) 
with E — h^ik^ — /c|)/2mf^^ Similarly, and fc_ for 
s_(x,fca;) are defined by Eq. ((T31) with E = h^{k1 + 
fc|)/2m. With suitable boundary conditions, we solve 
the differential equation Eq. H20|l to determine the semi- 
classical, one-electron spin densities. The total, spin- 
resolved, spin densities follow by inserting these one- 
electron quantities into 



s±(x) 



/±(k- 



cEt 



/±(k) 



s±{x:kx)-^d^\^ 



(21) 

This equation differs from Eq. ^ by the weighting fac- 
tor kx/ (k)^ This factor guarantees that the flux carried 
by each electron is proportional to its velocity (see Ap- 
pendix I). This is confirmed by Fig. 12b) which shows 
quantitative agreement between a fully quantum calcu- 
lation of s(x) using Eq. (fTU)) . Eq. ifT^ and Eq. © and a 
semi-classical calculation using Eq. H2()(l and Eq. H21(l . 

In light of the foregoing, it is reasonable to calculate the 
semi-classical single-electron spin current density from 



Q±{x,k,j;) = s±{x,kx) 



hkx 



m 



(22) 



and use the second equation in Eq. (j^ to find Q±(a;). 
The correctness of this prescription is illustrated in 
Fig.EIc). 
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FIG. 3: Cartoon of a spin spiral where the magnetization (ar- 
rows) rotates uniformly in the x-z plane of a fixed coordinate 
system where M(a;) ■ z = cos 6'. The inset shows a local coor- 
dinate system where M(3;) always points along the z axis. 



III. SPIN SPIRAL 

As a preliminary to our discussion of domain walls, it 
is instructive to discuss the spin density and spin current 
density for a spin spiral — an infinite magnetic structure 
where the direction of the magnetization rotates contin- 
uously as one moves along a fixed axis in space. Spin spi- 
rals occur in the ground state of some rare earth metalsSA 
and also for the 7 phase of ironiS^ Here, we focus on a spi- 
ral with uniform pitch p where the magnetization rotates 
in the x-z plane, i.e., in Eq. (j^J, 



and 



4>(x) 



0. 



(23) 



A cartoon version of this M(a;) is shown in Fig. O This 
figure also defines a local coordinate system that will be 
useful in what follows. The system (x',?/',z') rotates as 
a function of x so the magnetization M(a;) always points 
along -l-z'. 

CalvoS& solved Eq. (|10|l to find the eigenstates and 
eigen-energies of this spin spiral. In our notation, 



£±(k) 



and 



where 



and 



■!/'±(k, r) — e 
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(24) 

(25) 
(26) 

(27) 



From these results, it is easy to compute the single- 
electron spin densities defined in Eq. ifT^ . In the local 

(x', y' , z') frame. 



s'±{x, kx) = ±(0, sin a, cos a). 



(28) 



The corresponding calculation of the single-electron 
spin current densities Eq. (|22|l is straightforward but te- 
dious and not very illuminating. Therefore, we pass di- 
rectly to the total spin current density calculated by sum- 
ming over all electrons as indicated in the second line of 




FIG. 4: Electrons (wave vector k^) and holes (wave vector 
—kx) move in an effective field that is the sum of the ex- 
change field Bex (a:) || z' and a fictitious, velocity-dependent 
"gradient" field induced by the spatial dependence of the ex- 
change field. The spins align to the total effective field in an 
infinite spin spiral. The x axis lies in x'-z' plane. 



Eq. jnj. Again in the local {x',y',z') frame, 

Q'(x) =A(p,A:b)(0,0,1). (29) 

where A(p, fee) is a constant. This shows that Q(a;) oc 
M(a;), i.e., the spin current density for a free electron spin 
spiral is perfectly adiabatic. Wessely et alM- found con- 
sistent results in their density functional calculation of 
the steady-state spin current density associated with the 
helical spin density wave in erbium metal. We emphasize 
that Eq. H29|l is independent of pitch for an infinite spin 
spiral. As we discuss below, a similar independence does 
not hold for domain walls. In that case, wide walls are 
adiabatic, but narrow ones are not. 

The semi-classical formula Eq. ()22|l provides an appeal- 
ing way to understand the adiabaticity of the spin current 
density in the spin spiral defined by Eq. The key 

point is that the angle a in Eq. H26|l which fixes the direc- 
tion of s-|-(a;, kx) in Eq. (|28|l is positive when k^ is positive 
and negative when k^ is negative (Fig. ^J. Moreover, for 
every k^ electron that contributes to the shifted Fermi 
surface sums in Eq. @, there is a contribution from a 
—kx hole. Now, a hole has opposite spin density to an 
electron and the spin current density Eq. H22|) contains 
an additional factor of kx. Therefore, the two spin den- 
sity vectors in Fig. 0] subtract to give a net spin density 
along y while the corresponding two spin current den- 
sity vectors add to give a net spin current density along 
z'. This occurs for all kx in the sums so Q{x) aligns 
exactly with the local exchange field and thus with the 
local magnetization. 

The opposite situation occurs for fully occupied states 



4 



below the Fermi energy. The spins of the forward and 
backward moving electrons combine to produce a net mo- 
ment aligned with the exchange field, as necessary for 
self-consistency. Further, the spin currents, with the ad- 
ditional factor of kx add to give a net spin current along 
y, so that its gradient gives the correct form of the phe- 
nomcnological exchange torque density. 

To summarize: an electric current that passes through 
a spin spiral generates a spin accumulation with a compo- 
nent transverse to the magnetization. The spin current 
density possesses no such component due to pairwise can- 
cellation between forward and backward moving spins of 
the same type (majority or minority). Moreover, since 
the cancellation occurs within each band, the final result 
is insensitive to the details of intraband scattering. 

It remains only to understand the origin of the mis- 
alignment angle a. Why does each spin not simply align 
itself with Bex? Berger«4 was the first to notice this fact 
and the physics was made particularly clear by Aharonov 
and Sterni^ These authors studied the adiabatic ap- 
proximation for a classical magnetic moment that moves 
in a slowly varying field B(a;). Not obviously, the mo- 
ment behaves as if were subjected to a effective magnetic 
field Beff(a;) that is the sum of Bex(a;) and a fictitious, 
velocity-dependent, "gradient" field ^g{x) that points in 
the direction VB(a;) x B(x). For our problem. 
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(30) 



The presence of this gradient field is apparent from 
Eq. where the square root is proportional to |Boff|. 
The adiabatic solution corresponds to perfect alignment 
of the moment with Bcff(2:). This alignment is indi- 
cated in Eq. and in Fig. 0] More generally, the 
expected motion of the magnetic moment is precession 
around Boff (a;). Nevertheless, as indicated above, the to- 
tal spin current density for the spin spiral aligns with 
Bex(a;) (which is the conventional definition of adiabatic- 
ity for this quantity) when the net effect of all conduction 
electrons is taken into account. 



IV. DOMAIN WALLS 

Our main interest is the spin-transfer torque associated 
with domain walls that connect two regions of uniform 
and antiparallel magnetization. A realistic wall of this 
kind can be described by Eq. Q with22, 



9{x) — tt/2 — arcsin [tanh (x/w)] 



(31) 



The wall is Neel-type if 4>{x) — and Bloch-type if 
(j){x) = 7r/2. We will speak of the domain wall width 
w as "long" or "short" depending on whether w is large 
or small compared to the characteristic length 
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FIG. 5: Distributed spin-transfer torque for a long Neel do- 
main wall with w = 50 and L — 6.25 (kp = 1 and fcs = 0.4): 
semi-classical calculation of Nst (solid dots) compared to 
Eq. 13311 for Nad (a;) (solid curve). 



Intuitively, the adiabatic approximation should be 
vaHd when w ^ L. When applied to Eq. ((HJ, the pre- 
dicted adiabatic spin-transfer torque for our model is 

n TIP 

Nad(x) = -\ d^mix), (33) 

Z m 

where n is the electron density, and r] is the polariza- 
tion of the current. The calculations required to check 
this for long domain walls are difficult quantum mechan- 
ically (for numerical reasons) but straightforward semi- 
classically. At the single-electron level, adiabaticity again 
corresponds to alignment of the spin moment with the ef- 
fective field defined in Eq. (|3()|l . The results for a typical 
long wall (Fig. |3J) demonstrate that summation over all 
electrons produces alignment of Q(x) with M(x) so the 
adiabatic formula Eq. H33() is indeed correct in this limit. 

For short walls, we have carried out calculations of 
Nst (2:) both quantum mechanically and semi-classically. 
The two methods agree very well with one another (see 
Fig. 13) but not with the proposed form Eq. ^ . Bearing 
in mind that, when the magnetization changes, x' points 
along 9a;M and y points along M x SajM, our result for 
the spin-transfer torque is 

Nst(x) = Nad(a;) + a{x)S^ + b{x)y. (34) 

Nst (a;) differs from Nad (a;) because gradients in the 
gradient field induce single electron spin moments to 
precess around Bcff(a;) rather to align perfectly with it. 
Fig. |51 shows a{x) and h{x) as calculated for a typical 
short domain wall. The associated torques lie in the 
plane of the magnetization and perpendicular to that 
plane, respectively. These non-adiabatic contributions 
to the torque are both oscillatory functions of position 
that do not go immediately to zero when the the magne- 
tization becomes uniform. In other words, a{x) and b{x) 
are generically non-local functions of the magnetization 
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FIG. 6: Distributed spin-transfer torque for a short Neel do- 
main wall with w = A and L = 6.25 {kp ~ 1 and fc_g = 0.4): 
in-plane piece a{x) (heavy solid curve); out-of- plane piece b{x) 
(dashed curve); adiabatic prediction (light solid curve); sec- 
ond term in Eq. Q scaled to match the maximum of b(x) 
(light solid curve). 




FIG. 7: Non-adiabaticity in Eg. 1351 versus wall width scaled 
by the characteristic length in Eq. 13211 . Note the logarithmic 
scale. 



M(a;). The positive- valued function that falls to zero at 
the edges of the domain wall (light solid curve in Fig. ^ 
is the second function in Eq. ^ with C2 chosen to match 
b{x) at their common maximum. Evidently, the proposed 
torque function Eq. |^ gives at best a qualitative account 
of the out-of-plane non-adiabatic torque. 

A convenient measure of the degree of non-adiabaticity 
of the spin-transfer torque is 



Q 



max I b{x) \ 
max I Ni,d{x) \ 



(35) 



Fig. [7| plots this quantity as a function of scaled domain 
wall width w/L on a log scale. The observed exponential 
decrease of the non-adiabatic torque as the wall width 
increases can be understood from the work of Dugaev et 
alm^ These authors treat the gradient field in Eq. (|30|l as 
a perturbation and calculate the probability for an elec- 
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FIG. 8: Dependence of 7 in Eq. 1361 on domain wall sharp- 
ness: Squares are calculated points. Straight line is a guide 
to the eye that passes through the origin. 



tron in a (fc^; t) state to scatter into a (fc^ J,) state in the 
Born approximation. If we choose and k'^ as kp and 
kp, respectively, their results imply that the probability 
V that a majority electron retains its spin and becomes 
a minority electron as it passes through a domain wall is 



V cx exp(— 7w/L), 



(36) 



where 7 is a constant of order unity. This rationalizes 
the result plotted in Fig. because the magnitude of the 
minority spin component determines the amplitude of the 
spin precession around Boff(x) and thus the magnitude 
of the non-adiabatic component of s and Q in Eq. H22|l . 
In fact, A^ad oc 1/w, so it is the case that 



max I b{x) \ cx — exp{—^w/L). 
w 



(37) 



The slope of the straight line in Fig. [7| i.e., the value 
of the constant 7 in Eq. (|36|l depends on the sharpness of 
the domain wall. Using Eq. (|31|l and other simple domain 
wall profile functions, it is not difficult to convince oneself 
that a suitable measure of domain wall sharpness is the 
maximum value of the second derivative 6" (x) for walls 
with the same width. The numerical results shown in 
Fig. |S| confirm this to be true. The sharper the domain 
wall, the less rapidly the non-adiabatic torque disappears 
with increasing domain wall width. 



V. RELATION TO OTHER WORK 

1. Waintal & Viret 

Waintal and Viret^ (WV) used a free-electron Stoner 
model and the Landauer-Biittiker formalism to calculate 
the spin transfer torque associated with a Neel wall with 
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magnetization Eq. ^ and 

{0, a: < ~w 

{tt/2) [x/w + 1) , -w <x<w 
TT, X > W. 

(38) 

For this wall profile (which is exactly one half-turn of a 
uniform spin spiral in the interval ~w < x < w), WV 
reported oscillatory non-adiabatic contributions to the 
torque similar to our functions a{x) and h{x). This con- 
trasts with the perfect adiabaticity we found in Sec. IIIII 
for the infinite spin spiral. Moreover, the amplitude of 
the non-adiabatic torque reported by WV for this wall 
decreases only as \/w rather than (1/w) exp(— 7w/L) as 
we found above. 

The disparities between Ref. and the present work 
all arise from the unphysical nature of the domain wall 
Eq. Specifically, the divergence of 6"{x) at cc = ztw 

locates this wall at the origin of Fig. |H| where 7 = 0. This 
brings their result into agreement with Eq. (|37|l . Any 
rounding of the discontinuity in slope at x = ±w would 
yield a finite value for 6"{x) and thus a non-zero value of 
7- 

In Appendix II, we calculate the spin-transfer torque 
for the wall Eq. ^ using our methods. Qualitatively, 
the pure l/w behavior of the non-adiabatic torque comes 
from the fact that there is a sudden jump in 9'{x) at 
X ~ ibui. There is a corresponding jump in the direc- 
tion of BefT(a:) as defined by Eq. l|5n)l. Spins propagating 
along the the a;-axis cannot follow this abrupt jump and 
thus precess around the post-jump field direction with an 
amplitude determined by the sine of the angle between 
the before-and- after field directions. The latter is pro- 
portional to the jump in 9'{x), which is 'k/2'w for the 
waU Eq. (EHIl- 

2. Zhang & Li 

In spin spirals and long domain walls, we find that the 
non-equilibrium spin current is adiabatic, i.e., Q(a;) is 
aligned with M(a;) [or Bex(a;)]- At the same time, we 
find in both cases that the non-equilibrium spin density 
s{x) is not aligned with the magnetization; there is com- 
ponent of s(a:) transverse to M(a;). The corresponding 
transverse component of the spin current density can- 
cels between pairs of electrons moving in opposite di- 
rections. Zhang and Li^ found exactly the same form 
of non-equilibrium spin accumulation (called dni{x) by 
them) using a phenomenological theory. They proposed 
that this non-equilibrium spin density relaxes by spin- 
fiip scattering toward alignment with the magnetization. 
Such a relaxation would produce a non-adiabatic torque 
of the form given by the second term in Eq. The cor- 
rectness of this predicted non-adiabatic torque depends 
on the correctness of the assumed model for relaxation of 
transverse spin accumulation through spin flip scattering. 



Zhang and Li assume a form for the rate of spin flip 
scattering, Sm/rsU that has been used successfully as a 
phenomenological description of longitudinal spin relax- 
ation in systems with collinear magnetization. While it 
is plausible to extend this form, as they do, to describe 
transverse spin relaxation in non-coUinear systems, our 
calculations indicate that it is not likely to be correct. 
Our reasoning is simplest to appreciate for a spin spiral 
with small pitch p. In this limit, Eqs. H26|l and 1)28(1 show 
that the transverse component of the spin for every elec- 
tron eigenstate is proportional to its velocity. This means 
that the majority band electrons contribute a transverse 
spin accumulation and an electric current that are pro- 
portional to one other. The same is true, separately, for 
the minority band electrons. This conclusion is indepen- 
dent of the details of the electron distribution. There- 
fore, for a fixed total current, it is impossible to relax 
the transverse spin accumulation without changing the 
longitudinal polarization of the current. No such change 
occurs in the model in Ref. 0, casting doubt on the valid- 
ity of the form of the spin fiip scattering assumed there. 

Microscopic considerations also argue against this form 
of the relaxation. As we have emphasized, the nat- 
ural basis for an electron spin moving though a non- 
coUinear magnetization is not along the local exchange 
field Bex (a;), but rather along a local effective field 
Bcff(a;), which includes the corrections due to the gra- 
dient of the magnetization [see Eq. H30() ]. Any spin that 
deviates from parallel or antiparallel alignment with the 
effective field will precess around the effective field, and 
on average will point parallel or antiparallel. Thus, we 
expect that there is no tendency for electron spins moving 
in a non-uniform magnetization to align themselves with 
the local exchange field Box (a;) by spin-flip scattering (or 
any other mechanism). Rather, the adiabatic solution 
is precisely alignment of their spins with the local effec- 
tive field Boff(x). Without further microscopic justifica- 
tion, we believe that the phenomenological form of spin 
fiip scattering assumed in Ref. |^ should not be used in 
systems with non-collinear magnetizations. Hence, this 
analysis argues against the existence of the resulting con- 
tribution to the "non-adiabatic" torque from spin flip 
scattering. 



VI. SCATTERING 

We do not explicitly treat scattering in any of our cal- 
culations. However, the distribution function in Eq. J^l), 
a shifted Fermi distribution, is an approximate solution 
of the Boltzmann equation in certain limits. First, the 
electric flcld must be small enough that the transport is 
in the linear regime. Then, the appropriate limits are 
determined by three important length scales, the Fermi 
wavelength, the mean free path, and the characteristic 
length of the structure, either the pitch of the spin spiral 
or the width of the domain wall. In all cases, we con- 
sider the limit in which the Fermi wavelength is short 
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compared to the mean free path. This hmit allows the 
description of the states of the system in terms of the 
eigenstates of the system in the absence of scattering. 
Different limits apply to the cases of domain walls and of 
spin spirals because the distribution functions are inter- 
preted differently for these two structures. 

We use the Boltzmann equation in two different ways. 
When the mean free path is much longer than the char- 
acteristic size of the structure, the distribution function 
describes the occupancy of the eigenstates of the entire 
system. This distribution function is independent of the 
spatial coordinate and we refer to this approach as global. 
In the opposite limit, the distribution function is spatially 
varying and describes the occupancy of eigenstates of the 
local Hamiltonian, which includes the exchange field and 
the gradient field. We refer to this approach as local, as 
the distribution function can vary spatially. 

For spin spirals, the distribution functions are shifted 
Fermi functions of the eigen-energies of the spin spiral. 
In the limit that the pitch of the spiral is much shorter 
than the mean free path, the shifted distribution given 
in Eq. JBJ is a solution of the global Boltzmann equation 
in the relaxation time approximation. The distribution 
function also becomes a solution in the opposite limit, 
where the mean free path is much shorter than the pitch 
of the spiral. In this limit, the Boltzmann equation is 
considered locally rather than globally. At each point in 
space the states are subject to the local exchange field, 
and the local gradient field. The distribution function is 
defined for states that are locally eigenstates of the sum 
of the fields. The local distribution function is given by 
the adiabatic evolution in the rotating reference frames 
of the distribution function specified in Eq. ©. In the 
limit that the pitch of the spiral goes to infinity, this dis- 
tribution function locally solves the Boltzmann equation 
in the relaxation time approximation. Thus, for spin spi- 
rals, the distribution function given in Eq. is a solu- 
tion in the limits that the mean free path is much greater 
than or much less than the pitch. We speculate that the 
corrections in between these limits are small. 

Domain walls are not uniform in the way that spin 
spirals are, so the distribution functions need to be given 
a different interpretation. For these structures, the dis- 
tribution function is determined from the properties of 
the states in the leads. For example, in the Landauer- 
Biittiker approach to this problem,^ scattering is ignored 
in the domain wall itself and confined to the "leads" ad- 
jacent to it (these leads are assumed to be "wide" and 
function as electron reservoirs) . An applied voltage is as- 
sumed to raise the energy of electron states in one lead 
relative to the other. Thus, in a formula like Eq. ©, the 
distribution function is shifted in energy rather than in 
velocity. 

We also do not treat scattering within the domain wall 
explicitly, but we assume that the wall is bounded by long 
leads that are as "narrow" as the domain wall region and 
have resistances per unit length that are comparable to 
that of the domain wall region. Thus, the distribution of 



the states approaching the domain wall region is similar 
to the distribution of states in an extended wire, i.e., to 
that given by Eq. For domain walls in long wires, the 
distribution function for left going states is determined 
by the right lead and for right going states by the left 
lead. With this interpretation, the distribution given in 
Eq. © is a solution in the limit that the scattering in 
the domain wall is weak, that is, the domain wall is much 
narrower than the mean free path. 

The distribution in Eq. © is also a solution in the limit 
that the mean free path is much shorter than the domain 
wall width. Since the Fermi wave length is much shorter 
than the mean free path, it is much less than the domain 
wall width. In this case, quantum mechanical reflection is 
negligible and the quantum mechanical states are closely 
related to the semiclassical trajectories. With a similar 
interpretation of the distribution function as was made 
for the spin spirals in this limit, the same conclusion holds 
for the domain walls. 



VII. SUMMARY & CONCLUSION 

In this paper, we analyzed spin-transfer torque in sys- 
tems with continuously variable magnetization using pre- 
vious results of Calvo^S, for the eigenstates of an infinite 
spin spiral and of Aharonov and SterniSS for the classi- 
cal motion of a magnetic moment in an inhomogeneous 
magnetic field. Adiabatic motion of individual spins cor- 
responds to alignment of the spin moment not with the 
exchange field (magnetization) but with an effective field 
that is slightly tilted away form the exchange field by 
an amount that depends on the spatial gradient of the 
magnetization. Nevertheless, when summed over all con- 
duction electrons, the spin current density is parallel to 
the magnetization both for an infinite spin spiral and for 
domain walls that are long compared to a characteristic 
length L that depends on the exchange energy and the 
Fermi energy. 

Non-adiabatic corrections to the spin-transfer torque 
occur only for domain walls with widths w that are com- 
parable to or smaller than L. The non-adiabatic torque 
is oscillatory and non-local in space with an amplitude 
that decreases as exp{—^w/L). The constant 7 is 
largest for walls with the sharpest magnetization gradi- 
ents. This suggests that non-adiabatic torques may be 
important for spin textures like vortices where the mag- 
netization varies extremely rapidly. 

Using microscopic considerations, we have also argued 
that the role of the gradient field to tilt spins away 
from the exchange field casts serious doubt on a recent 
proposal by Zhang and LiS- that a non-negligible non- 
adiabatic contribution to the torque arises from relax- 
ation of the non-equilibrium spin accumulation to the 
magnetization vector by spin flip scattering. We con- 
clude that, if the second term in Eq. truly accounts 
for the systematics of current-driven domain wall motion, 
the physics that generates this term still remains to be 
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identified. 

Finally, we have carefully discussed the role of scat- 
tering in this problem with particular emphasis on the 
approximation used here to neglect scattering within the 
domain wall itself but to treat the adjacent ferromagnetic 
matter as bulk-like. We argue that this approximation 
is valid in limits that either include or bracket the most 
interesting experimental situations and therefore is likely 
to be generally useful. 
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Appendix I: Semi-Classical Weighting Factor 

The weighting factor k^/ {k) used in Eq. (|21|l brings the 
amplitude of the dynamic (transverse) part of the semi- 
classical, one-electron spin density into accord with the 
corresponding quantum mechanical amplitude. This can 
be seen from a simple model problem that we solve both 
quantum mechanically and semi-classically. Namely, a 
spin initially oriented along the -l-x direction propagates 
from x = — ootoa; = oo through a magnetization that 
changes abruptly from M(a;) = M(1,0, 0) for x < to 
M(a;) = M(0, 0, 1) for x > 0. For x < 0, the eigenstates 
are 



(39) 



^+(x) = 



(40) 



and for x > 0, the eigenstates are 

If we choose the incoming state as 

^(x) = Vf(2;), (41) 

the reflection and transmission amplitudes for spin flip 
irii,tii) and no spin flip (rf-f,if-f) are determined by 
matching the total wave function and its derivative at 
X = 0: 

V'f +'"TT(V'f )* +7'n(V-'j;)* = t^iip^ + tn'^i 

(42) 

/ciVf - ^■TT^T(V'f )* - ^n^i(V'r)* = ^TT^tV'i^ +iTifciV'f • 

(43) 

It is straightforward to confirm that these equations are 
solved by 



i-2 - ^2 



4V2fctfc| 



6fc|fc| -I- kf 



rn = 



2fcT(fci-feT) 
+ 6fc|fc| + kj 

2V2fct(fc| + k^) 
ki + 6fc|fc| + kl ' 



(44) 



0.2 



0.4 



0.6 



0.8 



kg /hp 



FIG. 9: The semi-classical weighting factor for the spin den- 
sity. Solid (dotted) curve is the expression on the left (right) 
side of the « symbol in Eq. 14811 . 



We are interested in the transmitted wave function. 



t^^l})Ux) + t||7/)+(x) 



which carries a spin density, 



(45) 



(x) = - [2int|i cos((5fcx), 2i||t|x sin((5fcx), [t]^ ~ t]{)\ , 

(46) 

where dk = k^ — k^. Notice that the oscillation is trans- 
verse to the X ^ oo magnetization and of amplitude 

If we analyze the same problem semi-classically, a ma- 
jority electron propagates freely until it reaches x = 0. At 
that point, the electron feels a magnetization perpendicu- 
lar to its magnetic moment and begins precession around 
that magnetization with unit amplitude. The associated 
spin density is 



SjJ:(x) — — [cos((5fcx), sin((5fcx), 0] 



(47) 



Comparing Eq. H46|l to Eq. (|47|l shows that the trans- 
verse oscillation amplitudes will be equal if we multiply 
the semi-classical result by the weighting factor 



32/c2fcx(fcT+A:i) 



2kf 



{k'^ + 6k^ki + kfy k^+ki (fc) 



= TTT- (48) 



Fig. 1^1 illustrates the quality of the approximation in 
Eq. H48|) if we identify fc| and k^ with kp and kp (re- 
spectively) in Eq. Q. Of course, k^ plays the role of k^ 
in Eq. (|?T|l. 



Appendix II: Spin Spiral Domain Wall 

The semi-classical spin density associated with electron 
propagation through a magnetization like Eq. ^ is the 
solution of Eq. (|20|l with suitable boundary conditions. 
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Choosing = 0, we simplify the notation by using the 
prefactor A = kg/{k) and an overdot for d/dx to write 
the components of Eq. (|2()|l as 



Sy 

Sz 



~ ~XSy COS 

= — Xsz sin 9 + Xsx cost 
— Xsy sinO. 



(49) 
(50) 
(51) 



In the local frame {x',y',z') defined in Fig. 13 the com- 
ponents of the spin density, 



Sx COS t 



Sx smO + Sz cost 



satisfy 



-As' - s'J 



Eliminating s' gives 



(A^ 



J n 



= 0. 



(52) 
(53) 
(54) 



(55) 
(56) 
(57) 



(58) 



The differential Eq. 1)58(1 cannot be solved analytically 
for realistic domain wall profiles. However, it is easily 



solvable for the wall defined by Eq. H38|l where one-half 
turn of a spin spiral with pitch p = tt /2w connects two 
regions with uniform (but reversed) magnetization. In 
the limit tt/ui <C A of a long wall, the components of the 
spin density transverse to the wall magnetization for the 
range x G [—w,w\ are (after multiplying the weighting 
factor kx/{k) for the semi-classical approach) 



S'xix) 
S'yix) 



2 (fc) 2wX 

h kx TT 

2 (fc) 2^ 



sin(A(a; ± w)) 

[1 — cos(A(a; ± w))] 



(59) 



where the plus (minus) refers to electrons that flow from 
left (right) to right (left). The associated spin current 
density and spin transfer torque carried by each electron 
follow from Eq. and Eq. {7)), respectively. Bearing in 
mind that x' varies with x, our flnal result for the torque 
(in the local frame) generated by a single electron moving 
from right to left is 



N' = 



N' 

y 



2~mjk}2w ^ 



cos A(a: — a)] x' 



' (fc) 2w 



sin A (a; 



a) y- 



(60) 



This may be compared with the results of Ref. |7i which 
pertain to the entire ensemble of conduction electrons. 
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